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Abstract 



We compute the dynamical structure factor S(q, r) of an elastic 
medium where force dipoles appear at random in space and in time, 

*73 ' due to 'micro-collapses' of the structure. Various regimes are found, 

depending on the wave vector q and the collapse time 9. In an early 
time regime, the logarithm of the structure factor behaves as (<7t) 3 ' 2 , 
as predicted in M using heuristic arguments. However, in an inter- 
mediate time regime we rather obtain a (<?t) 5 / 4 behaviour. Finally, 

rN | the asymptotic long time regime is found to behave as q 3 ' 2 r. We also 

give a plausible scenario for aging, in terms of a strain dependent en- 
ergy barrier for micro-collapses. The relaxation time is found to grow 

4/5 

with the age t w , quasi-exponentially at first, and then as tj, with 
logarithmic corrections. 



1 Introduction 

There is currently a great interest in soft glassy materials, such as pastes, 
foams or colloidal gels M. These systems are often out-of-equilibrium and 



exhibit interesting aging effects, typical of other glassy systems such as spin- 
glasses or polymer glasses || |J. For example, their rheological properties 
are found to depend strongly on the age of the system || |5|. From a the- 
oretical point of view, aging is expected both in response functions (such 
as magnetic susceptibilities or elastic moduli) but also in correlation func- 
tions, which are notoriously harder to study experimentally. Recently, a new 
experimental 'multi-speckle' technique has allowed one to properly investi- 
gate aging in a dynamical correlation function 0. The difficulty is to make 
the experimental acquisition time much smaller than the typical age of the 
system, otherwise the age of the system changes significantly during the ex- 
periment and the measured correlation function is meaningless. Using this 
technique, the dynamical structure factor of a colloidal gel made of aggre- 
gating polystyrene particles was measured and revealed several unexpected 
features |1J. (i) For a given waiting time t w (counted from the moment when 
the gel is formed), the dynamical structure factor S(q,r,t w ) is found to de- 
cay as exp(— A(qr) 3 / 2 ). This must be contrasted with the usual diffusive 
decay, as exp(— Dq 2 r): both the q and r dependence are anomalous, and 
non intuitive. The exponent 3/2 shows that the time decay is faster than 
an exponential, whereas one would have expected a slow decay in a glassy 
system. A similar exponent 3/2 was also found in different systems, such as 
diblock polystyrene/polyisoprene copolymers H, laponite || or other sys- 
tems || . (ii) The scaling in qr suggests that some kind of convection, rather 
than diffusion, is present in the system. However, since the system is den- 
sity matched, this cannot be a global sedimentation effect, (iii) When the 
age of the system increases, the relaxation of the dynamical structure factor 
becomes slower but retains its shape. In other words, only the coefficient A 
in the above expression is found to be age dependent. The corresponding 
relaxation time r r (q,t w ) (defined as A(qr r ) 3 ^ 2 = 1) is found to increase expo- 
nentially with t w at first, and then as tfa with fi < 1 for larger waiting times. 
Such an exponentially growing relaxation time with age was also reported 
for laponite |10| , but is very unusual. The second regime, however, is typical 



of most experimental glassy systems |TT], |T2"| |. 

A heuristic interpretation of the (gr) 3 ^scaling was proposed in |1|. The 
basic mechanism is that the gel contracts in a very heterogeneous way: lo- 
calized 'micro-collapses' create a long range elastic deformation field which 
is ultimately responsible for the g 3//2 behaviour. In order to get the correct 
r dependence the strain within these collapsing regions was postulated to be 



linear in r. This is however awkward since the instant r = has no special 
meaning and a linear growth would eventually lead to unbounded strains. 
Furthermore, this argument does not account for aging. 

The aim of this paper is to put the arguments of on a firmer footing. 
We show that an important time scale was left out in the analysis of Jj]], 
namely the typical collapse time which we will call 9. We find that for 
r <C 9, the relaxation of the structure factor for a given t w does indeed only 
depend on qr, but exhibits two distinct regimes: a short time regime where 
we recover the (gr) 3//2 behaviour, and an intermediate time regime where the 
decay is slower, as exp(— B(qr) 5 ^ 4: ), which might be of experimental relevance. 
For very large r's, the qr scaling breaks down and a new dynamical regime 
is found. We also give a simple argument to rationalize the dependence of r r 
with t,„. 



2 The model 

Following I], we assume that the dominant mechanism is the random ap- 
pearance of micro-collapses: since the micro-particles forming the gel attract 
each other rather strongly, the gel tends to restructure locally so as to create 
dense packing of particles. (Post mortem analysis indeed reveals that parti- 
cles actually tend to fuse together). Since the collapsing particles belong to 
a gel network, their motion will induce a certain strain field around them; 
other particles therefore move and dynamical light scattering probes this mo- 
tion. Assume for a moment that the gel is a one dimensional chain of beads 
and springs. When two particles decide to leave their equilibrium position to 
stick together, the left particle imposes a positive force +/o on the left part of 
the chain, while the right particle imposes an equal opposite force — /o on the 
right part of the chain. In other words, this creates a dipole of forces. This 
is also true in three dimensions, where the collapse of particles will result in 
a force dipole of intensity Pq in the direction n. When the dipole is formed 
at point r*o, the elastic strain field u at point r can be computed assuming a 
simple central force elasticity: 

KAu=-f(r), (1) 



where K is a compression modulus and f(f) is a local force dipole of the 
form: 

/(r) = /o[e-V5(f-f )]n, (2) 

where we will assume that e = en, i.e, the mean displacement of the parti- 
cles creates a force in the same direction. This dipolar force can be simply 
expressed by its Fourier transform: 

f k = iP (k.n) n exp(-ik.r ) (3) 

with P = f e. Note that a more refined model with shear modulus could 
be also be considered, but would only change some numerical factors in the 
following calculations. The solution of equation (HD is of course: 

u(f) = -— — — — n. (4) 

AttK \r -r | 3 

The r~ 2 dependence of the strain field has an immediate consequence which 
will be of importance in the following: if there is a finite density of force 
dipoles randomly scattered in space, the probability that the stress at a given 
point has an amplitude u decays for large u as w~ 5 / 2 , which has a diverging 
variance. (This divergence is however cut-off if the finite size of the dipoles 
is taken into account). This property of the distribution of displacements 
and stresses will be responsible for the unusual g-dependence of the structure 
factor. Note also that from dimensional arguments, Pq/AttK = v o is a volume 
given by £ 3 A 2 , where £ is the typical size of the collapsing region, and A < 1 is 
the contraction ratio (that we will assume to be of order 1 in the following) . 
Now, let us assume that the micro-collapses are not instantaneous but 
take place progressively, over a time scale 9. A given event j starts at time tj 
and is completed at time tj + 9; the dipole intensity P(t) at time t is a certain 
function of (t — tj)/8. We will assume that this function is approximately 
linear, and simply write Pj{t) = Po(t — tj)/6, which saturates at Po when 
the collapse is completed. (Note that a micro-collapse could itself be the 
result of many successive events). The dynamics of the individual particles 
is presumably dominated by viscous friction, therefore we write the following 
equation of motion for the strain field u: 

1 ^&r = K ^ + T. m t) + v(r, t), (5) 



where the Fourier transform of the dipolar force is 

fj(k,t) = iPjit) (k.Hj)nj exp(-ik.fj) (6) 

and 7 is a friction coefficient, rf is the thermal random force, and we now take 
into account the fact that many microcollapses take place, at different times 
tj and different positions fj, with different orientations Hj of the dipoles. 
Actually, we will assume in the following that these events occur randomly 
in space and time, with a certain rate p per unit volume and unit time (see 
below for a further discussion of this assumption). The quantity K/j is a 
diffusion constant that we will call D. Equation (|J) defines the model that 
we want to study and from which we will compute the dynamical structure 
factor S(q,r), defined as: 

S(q, t) = (exp [iq ■ (u(f, t + r)- u{f, £))]) (7) 

where the brackets refer to a spatial average over r or, equivalently, over the 
random location and time of the micro-collapse events. In the following, we 
will neglect the thermal random force, which would add a Debye- Waller diffu- 
sive contribution to the dynamical structure factor, and set rf = 0. However, 
as we discuss below, the presence of rf has an indirect crucial effect since the 
nucleation of micro-collapses is probably of thermal origin. 

3 The slow collapse regime 

A first step is to calculate the Fourier transform of the time derivative of the 
displacement field u(f, t) created by a single dipole located at fj, in direction 
Hj, that we denote v(k,t\fj,fij,tj). One finds: 

v{k,t\rj,tj) = -rexp(-rk ■r j )—— L J exp(-Dk t) 

\exY>{Dk 2 tj) - exp(Dk 2 min(i i + 9, £))] . (8) 

The displacement difference between t and t + r can therefore be expressed 



as: 



'•*+ T , , ^ r d 3 k 



u(r, t + r) - (r,t) = J dt' Y. J 77p^ exp(ik ■ r)v(k, t\fj, tj) (9) 



i/tj<t 



The analysis of this expression reveals that there are a priori five different 
cases to consider for the relative position of the time tj when the j th micro- 
collapse takes place and the other relevant times: (al) tj < t — 9 < t' — 9: 
the j th event is over before t, t', (a2) tj < t' — 9 and t — 9 < tj < t + r — 9; 
and (bl) tj > if - 9, t - 9 < tj < t and tj < t + r - 9, (b2) tj > if - 9 and 
t + r -9 <tj <t and finally (b3) tj > tf - 9 and t < tj < t + r. 

We first consider the 'slow' case where the experimental time r is very 
small compared to the collapse time 9. The contribution of the different 
regimes to S(q, r) can be estimated and one finds that a new, g-dependent 
time scale r q appears, defined as: 

DO „ 



Tn 



qv 



(10) 



such that, depending on the ratio T/r q , the dominant contribution to the 
decay of S(q,r) comes from different regions of the t',tj plane. Let us first 
consider the case r -C r q . The dominant contribution then comes from region 
(b2). Introducing Tj = t — tj, the average contribution of a event i to S(q, r) 
is given by: 



/t dv f f dn I 

— — / d(cos a) / — < exp 



i[q-n)cosa— —J r {- 



9r 2 2, AD 



- 1 



(11) 
where V is the total volume of the sample, and T is a function defined as: 

T{x) = erfc(x) - x erfc'(x) J"(0) = 1. (12) 



Now, introducing in Eq. (|Tl|) r = sJ\q ■ n\TV /9 and Tj 
one can rewrite (JTT]) for very large volumes V as: 

3/2 

dn I / — — <; =- sin 



exp 
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with R = T\q-n\/r q z. Each micro-collapse contributes independently; and 
since within a small time interval 9dz one has a total of pV9dz events, the 
contribution of the (b2) regime reads: 

3/2 

/ - - / . _. ; i'/ • " i / in i 
exp 




Since we have assumed that r <C r q , it is justified to set R = provided all 
integrals (over s and z) converge. One can check that this is indeed the case. 
The final result reads (for t > 9): 



S(q, r) = exp 



75~ 



p6{D6) 



3/2 



3/2' 



(t«t,), (15) 



which has the form suggested by the arguments of Q, in particular, it indeed 
only depends on (gr) 3 / 2 . (The numerical value of the prefactor is 1.67996...) 
However, the dimensional factors differ from those obtained in Jl|, and the 
result is only valid in the short time regime r ^ 9 and r <C r q . Note that 
the combination p = p9(D9) 3 ' 2 is adimensional and represents the average 
number of events taking place within a time interval 9 and within a diffusion 
volume (D9) 3 ' 2 . The contribution of the other regions in the t',tj plane can 
be analyzed similarly, and lead to sub-dominant corrections proportional to 
P( r / r <?) 2 an d f>( T / T q) 5 ^ 2 - Therefore, if rjr q is not very small, an effective 
exponent somewhat larger than 3/2 can be observed. 

Turning to the regime r 9 < r <C 0, we now find that the dominant regime 
is (al). Introducing Tj = t — 9 — tj > 0, one self-consistently finds that the 
relevant region is Tj *^> 9. This enables one to write the average contribution 
of a single event as: 



exp 



D9 2 



T 



3/2 



T Q T J 



dcosal cosal 3 ' 2 



s 2 ds 



v {g^vJv) 



sin 



G^VRf) 



with now cos a = q ■ n/q, Rl = t9 2 \ cos a\/r q r 2 and: 

2 



G{x) 



3 -x 2 

x e . 



The adequate change of variables is now found to be Tj = z9Jr/r q , which 
means that Tj is indeed much larger than 9 for most values of z. All the 
integrals (over z, n and s) can be performed exactly in the limit r q <r and 
the final result reads: 



S(q, t) = exp 



2 10 7r 3/8 r( _| )r( ^ )sinf 
34517/8 



T 



Tn 



5/4' 



[T q ^T < 



(18) 



where the ridiculously complicated numerical factor is equal to 1.00993... 
Therefore, we find that in this regime the qr scaling still holds, but the 
power 3/2 is replaced by 5/4. Again, the other regions of the plane t', tj lead 
to sub-dominant corrections, the most important being oc p(r/Tg). 

Before discussing the 'fast' collapse regime, let us insist on one further 
limitation of the above calculation: we have assumed that the dipolar field 
behaves as 1/r 2 everywhere in space, that is, that the dipoles are point-like. 
This is valid is q is large enough. If q is very small, the dominant source of 
the decay of S(q, r) comes from points were the displacement is maximum, 
i.e. in the immediate vicinity or 'inside' the collapsed regions, where the 
whole analysis breaks down. 

4 The fast collapse regime 

Let us now study the regime where r ^> 6, corresponding to 'fast' micro- 
collapses. (Note that the fact that 6 is small does not mean that the events 
are frequent: this is described the nucleation parameter p). The calculations 
are very similar to the above case. The relevant time scale which now appears 
naturally is 

f q = qv /D = 6 2 /r q . (19) 

When t <C f q , we find exactly the same decay as Eq. (|T8| ) above: note 
indeed that 6 actually drops out of this expression. However, when r ^> f q , 
the dynamical structure factor reads: 



S(q, t) = exp 



16V2vr3 3/2 , 



75 



(r»f,), (20) 



which is independent of the diffusion constant D. Therefore, the asymp- 
totic decay of S(q, r) is a pure exponential, with an anomalous decay time 
oc g -3 / 2 . This is the result one obtains if all retardation effects are neglected 
(i.e. D — > oo): the factor pr simply counts the average number of events 
between t and t + r, and g 3//2 reflects the fact that the distribution of local 
displacements u decays as w~ 5 / 2 and has a diverging variance (see, e.g. fll3| ). 
For a distribution with a finite variance, one would obtain the usual q 2 de- 
pendence. The appearance of power-law tails in systems with long-ranged 
interactions is well known. It is called the Holtsmark distribution in the 



context of gravitational fields (T^| 



5 Summary of the results and discussion of 
the experiments 

It is useful to summarize our results in a schematic way, in terms of the 
behaviour of — log S(q, r). The two physical cases depend on the relative 
position of r q and f q , or equivalently on the ratio D9/qv . For DO <C qv one 
finds r q <C 9 <C f q and: 

(grp (r^r q ); for) 6 ' 4 (r g «r«f g ); q^r (r»f g ), (21) 

whereas for D# 3> gt>o, the (gr) 5//4 regime is squeezed out and the results are 
simply: 

{qrfl 2 (r<0); g 3/2 r (r>0). (22) 

It is easy to check that the above expressions (with their corresponding di- 
mensional prefactors) correctly match one another in the crossover regions. 

To which of these regimes does the experiment correspond to? First, a 
rather good qr scaling across most of the time regime (r < 10 4 seconds, 
or 3 hours). This means that the asymptotic regime is not observed, and 
therefore that max(#, f q ) > 10 4 s. Let us first suppose that the experiments 
are in the regime D9 -C qvo where r q <C 9 <C f q . 

According to our calculations, a (gr) 5 ' 4 behaviour should be observed for 
r q <^ r <^ f q . Fig. 1 shows that such a possibility is indeed compatible with 
the experimental data, although the data is by no means compelling. Under 
this assumption, the crossover time indicated in Fig. 1 should be interpreted 
as qr q = D9 2 /vq, which is thus of order 3. 10 7 cm^ 1 s. It is reasonable to 
assume that the volume of the collapsing region v o ~ £ 3 is of the order of 
the cluster size at gelation determined in 0, i.e. 3. 10~ 8 cm 3 (£ ~ 30/im). 
The elastic diffusion constant can be estimated from the elastic modulus of 
the structure [G' ~ 1CT 3 dynes/cm 2 ), the value of £ and of the viscosity 
i] ~ 10~ 2 Po: D = K/'y ~ G'v /r]^ ~ 10 -6 cm 2 /s. Therefore one gets 9 on 
the order of 1000 seconds. For the maximum experimental value of q ~ 6000 
cm -1 , one finds r q ~ 5000 seconds, which is a factor five larger than 9, in 
contradiction with our hypothesis. However this factor five is marginal in 
view of the roughness of our estimates. For example, the collapsing volume 
vq could be somewhat larger and numerical factors could help. 

Another possibility is that the (qr) 5 ' 4 regime is in fact not seen experi- 
mentally, i.e. D9 3> qvQ. This however requires that 9 is actually quite large, 

9 



at least 10 4 seconds. It would be interesting to determine p or 9 from an 
independent experiment in order to clarify this question. 

6 The aging phenomenon 

We finally discuss an important aspect of the experiments that we left out up 
to now, namely the fact that the dynamical structure factor is age dependent: 
the relaxation time grows with the age t w of the sample. This growth is even 
unusually (exponentially) fast in the initial stages fl], [K| . 

The basic mechanism for aging is that micro-collapses induce a tensile 
strain on the non collapsed regions. Therefore, these regions find it harder 
to collapse subsequently. We assume that the volume density of dipoles is 
4>(t w ). Therefore £* = -1 / 3 is the typical distance between two dipoles. Let 
us assume also that micro-collapse events are thermally activated, so that 
their rate p{t w ) can be expressed as: 

/ AE(t w ) \ 
p(t w ) = po exp I -j-jr- I , (23) 

where AE(t w ) is an energy barrier and T the temperature. The above mecha- 
nism means that the barrier height depends on the typical local strain | Vu\t yp . 
More precisely one can write that AE = y(^o + \^u\ty P £o} 2 — y^o, where 
£o is the equilibrium length of one elementary elastic unit, supposed here 
to be a constant. For small enough |Vw|ty P , one expects a linear relation- 
ship: AE ~ K£l\Vu\t ypi whereas for large | Vif | i yp this dependence becomes 
quadratic: AE ~ y^ol Vwl^. Obviously, new events will preferentially occur 
in less strained regions, i.e., far from collapsed regions, inducing some spatial 
correlations between micro-collapses that we have neglected above; hence, we 
will assume in the following that the small strain regime AE ~ K£q | Vm | t yp is 
the dominant one. The typical strain can be estimated as the sum of strains 
induced by all the dipoles surrounding a given particle in the system. For a 
unique dipole of orientation n at a distance r of the particle, one obtains: 

|V"| 2 = jt(3^ + U. (24) 

From this expression, we obtain | Vu|L„ = (|Vm| 2 ), where the brackets denote 
the average over all contributions of dipoles at a distance r > £* from the 

10 
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Figure 1: Structure factor S(q,r) as a function of qr, averaged over different 
q (on the order of 10 3 cm -1 ). The data here is the same as that of Fig. 2, 
(U, and corresponds to a waiting time of t w = 67 hours. Main figure: we show 
log(— log S(q, t)) as a function of loggr. The data is linear with a slope 3/2 for 
most of the time interval, but bends down at long times. This effect is more 
clearly seen in the inset, where we have plotted log S(q, r) as a function of (qr) 3 ' 2 . 



A (qr) 5 ' 4 behaviour is shown for comparison, and suggests that qr q 
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particle, the dipoles being distributed with a volume density = l/£* . The 
result is 

\Vu\typ = ( -jT- J V (j) (25) 

Hence, AE/ksT = T(f>, where T ~ p^fo- 

We want now to evaluate the age dependence of p{t w ) and 4>(t w ). The 
time evolution of <j) is given by: 

d(j> ( AE{4>) \ 

— = p(0)=p o exp^^-^j (26) 

The solutions can be written explicitly and one finds that the density of 
dipoles increases logarithmically with time: 

<KU = 0(0) + \ log (l + ^) , (27) 

where to is the initial activation time to — [p(t = 0)r] _1 . The rate of forma- 
tion of dipoles (micro-collapses) therefore decreases with time as: 

*•) = fs^) < 28 > 

The above calculations for S(q, r) can be simply extended to the case 
where p is not constant provided p does not vary too much on the scale 
of r, which is true if t w is sufficiently large. In order to compare with the 
experimentally determined relaxation time r r , one must solve the equation 
S(q,r r ,t w ) = e" 1 . Assuming that the experiments are in the (gr) 5//4 regime, 
one finds: 

r r oc (t w + to)" (29) 

with the exponent p = |. Note that the fact that p < 1 here is a conse- 
quence of the functional form of the non exponential relaxation. The basic 
mechanism, i.e. the nucleation a micro-collapses, follows a 'simple' aging law 
p ~ t" 1 . Correspondingly, if the experiments are in the (gr) 3 / 2 regime or 
in the final predicted regime in g 3//2 r, p would respectively have the values 
p = \ and p = 1. The best agreement with the experimentally determined 



r r ~ t°-» M is for p - i 
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Thus, as observed experimentally, r r seems to saturate to a finite value t$ 
for t w — > and to grow as t 4 / 5 for long times (see 0, Fig. 5). All this results 
have been obtained in the regime of small strains. If the strains become 
large (i.e. at large times), the energy becomes quadratic in |Vii| iOT and reads 

AE/ksT = T'(fi 2 , where V ~ k~~r v o- ^ n ^ ms case > logarithmic corrections 
appear in the previous expressions and one obtains in the large time regime: 
r r oc tf^ log M ' 2 [p V / T't t0 ]. More generally, if AE oc (j)", one finds an initial 
regime where r r grows as exp(t^) at first, and then as t^, with logarithmic 
corrections. 

Finally, the above mechanism suggests that when the density of dipoles 
becomes of order unity, the system should macroscopically collapse, as indeed 
seen experimentally fij at long times. 

7 Conclusion 

In conclusion, we have computed the dynamical structure factor S(q, r) of an 
elastic medium where force dipoles appear at random in space and time, as 
'micro-collapses' appear in the structure. The balance between these collapse 
events and the elastic relaxation of the internal stresses in the medium lead to 
various regimes, depending on the wave vector q and the collapse time 9. In 
an early time regime, the logarithm of the structure factor behaves as (gr) 3 / 2 
plus subleading corrections, as anticipated in using heuristic arguments. 
However, in an intermediate time that might be relevant for the experiments 
of |l|, we obtain a (gr) 5 / 4 behaviour that would have been difficult to guess 
from simple arguments. Finally, the asymptotic long time regime is found 
to behave as g 3 / 2 r, where the qr scaling is not obeyed; but this last regime 
seems not to be observed experimentally and it is likely that before that, a 
macroscopic collapse occurs as the one observed in colloidal gels 0. 

The relevance of the microscopic mecanism described in this paper still 
remains to be checked experimentally, in particular in order to know if the 
micro-collapses, if they indeed exist, continuously occur in time and space or 
if they only occur in a short time interval t c after the formation of the gel. 
We have made a calculation similar to the one above for this last case. In the 
limit where 9 <C t w , we find — log S(q,t w ,t w + r) oc t c (qr/t^ 4: ) 2 for r <^ t w 
which then saturates to a finite value (which goes to 1 as t w — > oo, since in 
this limit nothing moves any longer). 
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We have also given a plausible scenario for aging, in terms of a strain 
dependent energy barrier for micro-collapses. The relaxation time is found to 
grow with the age t w , exponentially at first, and then as tf^ 5 with logarithmic 
corrections, which seems to be in good agreement with experiments. 

It would be interesting to analyze other experiments where similar effects 
have been reported 0, @, 0, in particular in micellar crystalline phases and 
onion phases, along the lines of the present work. In these systems however, 
motion of defects such as dislocations or grain boundaries may also play a 
crucial role. 
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